Polarization of interacting bosons with spin 
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We prove that in the absence of explicit spin-dependent forces one of the ground states of in- 
teracting bosons with spin is always fully polarized. Generally, this state is degenerate with other 
states, but one can specify the exact degeneracy. For T > the magnetization and zero-field 
susceptibility exceed that of a pure paramagnet. The results are relevant to experimental work 
on triplet superconductivity and condensation of atoms with spin. They eliminate the possibility, 
raised in some theoretical speculations, that the ground state or positive temperature state might 
be antiferromagnetic. 

PACS numbers: 05.30.Jp, 03.75.Fi, 74.20.Rp 



Bosons with spin first appeared in condensed matter 
physics in the theory of 3 He superfluidity, which is due to 
spin-triplet p-wave pairing of the Helium atoms jj| . More 
recently they were discussed in connection with the rising 
interest in triplet superconductivity ^ |[ 0, || and mul- 
ticomponent Bose-Einstein condensation || [?], The 
renewed fascination with spin-triplet superconductivity 
was generated by recent experimental evidence for triplet 
pairing in heavy fermion systems j^j , organic conductors 
j|, and the layered oxide Sr 2 Ru04 j|, Recent ad- 
vances in the study of the Bose-Einstein condensation 
phenomenon enabled condensation of bosons with non- 
zero spin (spin-1 in particular), such as 23 Na and 87 Rb 
atoms, by confining them in optical traps rather than the 
usual magnetic traps, thus preserving their spinor nature 

§0- 

Moreover, it is possible to have condensates of bosons 
with other internal degrees of freedom which behave like 
spin, thus realizing "spin-1/2" bosons. For example, a 
two component Bose gas was produced in a magnetically 
trapped 87 Rb condensate by rotating two hyperfine states 
into each other, creating an ST/ (2) symmetry gj. 

All these experimental achievements have generated 
new interest in the study of interacting spinor bosons. 
In particular, some tantalizing hints for a possible con- 
nection between the onset of triplet superconductivity 
and the appearance of internal magnetic moments Q 
call for better understanding of the magnetic properties 
of strongly interacting bosons. We report here an exact 
and general result, stating that for a generic Hamilto- 
nian of bosons with spin, a fully polarized state is among 
the ground states - as long as there are no explicit spin- 
dependent forces - however complicated the many-body 
interaction potential might be. We can also classify all 
other ground states. The conclusion is that the ground 
state susceptibility is infinite. We further show that the 
positive temperature magnetization and zero-field sus- 
ceptibility exceed those of non-interacting distinguish- 
able (Maxwell-Boltzmann) spins. 

Recent works H analyzed the two body interaction for 



spin-1 atoms in an optical trap. Using an effective low 
energy approach, the interaction term was parameterized 
there by 



V(n -rj) = S(Ti - r 3 )(c + c 2 S 1 ■ S 2 



(1) 



where the coefficients cq and c 2 are related to the s-wave 
scattering lengths in the singlet and triplet channels. The 
possibility of an antiferromagnetic coupling c 2 > has 
attracted much attention in the past few years. Ex- 
perimental estimates for these scattering lengths suggest 
that c 2 > for a condensate of 23 Na atoms 10 . In 
such a case, the energy is minimized by a polar state 
with (S) = 0. We point out that in view of the follow- 
ing exact results this scenario can be ruled out as long 
as the full many-body interaction is spin- independent. 
In fact, as we explain later, the theorem below shows 
that the effective two-body spin-dependent interaction 
described in Eq. ([!]) can not describe the low-energy be- 
havior of bosons with spin-independent interaction, no 
matter what the sign of c 2 is. Therefore, a bosonic spin- 
spin term of the type (0) can only be justified by go- 
ing beyond a spin-independent bosonic Hamiltonian, ei- 
ther by introducing explicit spin-dependent forces, or by 
taking into account the underlying fermionic physics of 
the atoms. This situation should be contrasted with the 
fermionic case in which effective exchange interactions, 
ferromagnetic or antiferromagnetic, do come about from 
spin-independent Hamiltonians. 

The Hamiltonian of a Bose gas of N bosons, and with- 
out spin-dependent forces, is given by 



H 



N 1 

^ 2m 1 



nxi,x 2 , ■ 



,xat) 



(2) 



where m is the mass of the bosons, and xi,x 2 , • • • , Xjv 
are the spatial coordinates of the bosons, v is a spin- 
independent, totally symmetric, potential, which in- 
cludes any confinement, disorder, and A-body interac- 
tions. We look for the lowest energy state that is totally 
symmetric under permutations of the space-spin indices. 
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To find it, we temporarily ignore both the spin of the par- 
ticles and the identity of the bosons at first, and look for 
the ground state of N non-identical particles, i.e., parti- 
cles without spin and without any symmetry restriction 
on the wave-function. We seek the absolutely lowest state 
of H, V(xi,x 2 , • • • ,xjv). 

It is known (e.g., p|) that this ground state is unique 
(up to an overall phase) and satisfies 



V>(xi, • • - ,Xiv) > 



for all 



Xi,---,Xjv 



(3) 



(unless v is so repulsive that it prevents particles from 
changing places by a continuous path - which would be 
non-physical). This result holds for all boundary condi- 
tions (Dirichlet, von-Neumann, periodic). Eq. (H|) follows 
immediately from the fact that ip can be chosen to be real 
(because ip* is also a ground state and hence ip±ip* are 
ground states), and from the fact that the variational en- 
ergy of the function l^l then equals that of ip itself 
Eq. (6.17.2)]. Hence \i\)\ must be a ground state. From 
this we conclude that if ip is a real ground state, either 
tp = a\ip\, and thus (||) is satisfied, or \ip\ — ip is a ground 
state that is zero on a set of positive measure, which is 
impossible unless v has the pathology mentioned above. 
If, now, <j> is another real ground state that is orthogonal 
to ip then 4> — \4>\ up to a phase (as we just proved). 
But two positive functions cannot be orthogonal, so ip is 
unique. 

Equation (g) implies that ip must be totally symmet- 
ric. I.e., if we sum tp over all permutations we obtain 
a function that is (a) symmetric, (b) non-zero, and (c) 
a ground state. By uniqueness, this symmetric function 
must equal ip up to a constant. 

We now return to our identical spin-full boson problem, 
and define 



(xi, • • • , xjy; oi, • • • 



,crjv) = 1p(xi,- ■ - ,Xjv)| TTT ' ' •> 



This function is totally symmetric in both spin and spa- 
tial indices, and thus is a valid wave-function for bosons 
(with spin). As we have shown it minimizes the energy 
regardless of symmetry restriction, and is, therefore, a 
ground state of the boson system. 

This ground state has spin angular momentum 
J = NS, where 5 is each boson's spin, and carries the 
usual (25 N + l)-fold degeneracy. There will be other 
ground states, but each must be a spin function times the 
unique spatial function we have just described. The ques- 
tion of the ground state degeneracy thus reduces to the 
following question: How many different spin functions 
are there that are totally symmetric (up to the trivial 
2 J + 1-fold degeneracy)? The answer is as follows. 

(i) . If 5 = 1/2 ("spin-1/2" bosons @) then there is just 
one function, J = N/2, and the ground state is therefore 
non-degenerate (except for the trivial degeneracy). 

(ii) . If S = 1 there are functions with J = N, N — 2, 
N — 4, ... and each of these appears exactly once. Wave- 
functions with J = N—l, N—3, N — 5, ... do not appear 



in the ground state. This is contrary to what one would 
have for a paramagnet in which all J values are degener- 
ate, and J-values smaller than N occur multiple times, 
(iii). Other values of 5 can be similarly analyzed by 
studying representations of SU(2S + 1) or by studying 
(2S + l)-rowed Young's tableaux. 

If a magnetic field term Hm = — M^SiLi "! (where 
h is the field, and \io~\ is the atomic magnetic moment 
in the z direction, including the g- factor), is added to H 
then the ground state energy shifts by —/ihNS, i.e., the 
zero-field susceptibility is infinite. 

We now discuss the positive temperature case and com- 
pare it with pure paramagnetism. In the following, a 
pure paramagnet means a set of N non-interacting, dis- 
tinguishable, spin-5 particles, which has no degrees of 
freedom besides the spin. The pure paramagnet magne- 
tization is, therefore, 



Af para (T, h)= Nn 
_ Nji 



J2<r=-S (XeX -P[/ 3 ( J > h(X ] 



EL-s exp[/3/z/io-] 



(4) 



(25 + 1) coth 



Pfih 



coifi ( —0(ih 



where (3 = l/k B T. 

Our theorem about enhanced magnetization in the 
ground state can be generalized to positive temperature 
T, as follows. 

THEOREM: For each T and magnetic field h > 0, the 
magnetization M(T, h) is greater than the pure param- 
agnetic value M para (T, h). Moreover, the zero- field sus- 
ceptibility dM(T,h) I 'dh\h=o, also exceeds that of a pure 
paramagnet, dM pma (T,h)/dh\ h=0 . 

Proof. — Let us denote a particle configuration 
X\,X2, ...,xm by X and a spin configuration (in the o~ z 
basis) <7i, CT2, ...er n by E. Since spin does not enter H 
except through Hm the (operator) Boltzmann factor has 
the kernel 

K{X, X'; T) exp [/3/z/i5(E) ] 

where 5(E) = J2iLi °fi <fe,£' is the 'Kronecker delta' 
n^Li^.o-' and K is the (unsymmetrized) kernel of 
exp[— (3H ]. The symmetrized Boltzmann factor is then 



Sx,nK{X,X';T)exp[pfihS{T,)} fc )B , 



(5) 



and 5x,£ is the symmctrizcr on the X, E variables. (It is 
not necessary to symmetrize on the X' ,YI since this will 
be automatic.) 

The important point to notice is the T > analogue 
of (g), namely, K{X,X';T) > for all X,X' (outside 
of a 'hard-core' region). This can easily be seen by the 
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path-space (Wiener) integral for exp[—/3H](X,X') via 
the Feynman-Kac formula - or else via the Trotter prod- 
uct formula Q. 

The Trotter formula states that exp[ — /3H ](X, X') is 
the limit as N — > oo of 

J A(X,X 1 )B(X 1 ,X 2 )A(X 2 ,X 3 )B(X 3 ,X 4 )--- 

A(X 2 n-2, X 2 n-i)B(X 2 n-i, X')dXi ■ ■ ■ dX 2N _i (6) 

where A(X,Y) = exp[-/3H /N]{X,Y) (with 
Hq = kinetic energy operator) and B(X,Y) = 
exp[-0v/N](X,Y) = exp[-0v(X)/N]5(X -Y). The 
kernel A is a Gaussian, A(X, Y) = C\ exp[C 2 (X - Y) 2 ]. 
Since all kernels are nonnegative, we see that the multi- 
ple integral ((^) yields a nonnegative function. The limit 
of nonnegative functions is nonnegative and, with some 
work one can show that K(X,X';T) > for all X,X' 
(outside of a 'hard-core' region). In case exp[— v(X)j 
is not integrable we must first approximate it by a 
bounded function and then remove the approximations 
at the end. The positivity of the kernel is conserved 
under symmetrization. 

In fact, the only thing we really need is the weaker 
assertion that K V {X;T), defined below, is nonnegative 
and is positive on a set of positive X measure. 

To calculate the partition function Z(T, h) we set X = 
X' and £ = £' (after applying Sxx) an d then integrate 

over X and sum over £ (i.e., J^t, = E*Li Ect 1= -s)- We 
obtain an expression of the form 

Z ( T ' ft ) = / dXK v (X- : T) 

exp[/?^/i5(S)](5 7r(s)iS , 

s 

where the first summation is over all permutations tt and 
K 7T (X;T) is a 7r-dependent function of X and T defined 
by K„(X;T) = K(irX, X;T). 
To compute M(T, h) let us define 



Z„(T, h) = ^exp [(3(ihS(£)] 



M T (T, h) 



Es g(g) CX P [P^ h S{^)] 

Z„{T, h) 



C W (T, h) = ZAT,h) f K n (X;T)dX 



Then we can write 



M(T, h) 



Y^OaX h)M«(T,h) 
En CAT, h) 



(7) 



(8) 



(9) 



(10) 



min^ M n (T, h) . (The strict inequality > follows from the 
facts that C„(T,h) > and that the numbers M^(T,h) 
are not all equal.) 

Our problem reduces to deciding which tt will make (^) 
as small as possible. It is the identity permutation, and 
this gives M para (T, h) as the lower bound. To see this we 
note that (T, h) is the magnetization of distinguish- 
able particles that interact only through the constraint 
given by S w ^),^- If ^ is the identity, there is no con- 
straint, and one gets the pure paramagnetic magnetiza- 
tion (@). For any other permutation tt, ^(sj.s connects a 
group (or groups) of, say, k > 1 spins. There is one group 
of size k for each cycle of length k in tt. This grouping 
only raises the magnetization by reducing the entropy of 
the low moment states. In other words, the assertion 
follows from the inequality 

, Eff=-s0-exp[^u/ia-] EL-sNexp[M(m 

K T, < 



whose validity can be immediately seen by rewriting it 
as kM(T, h) < kM(T/k, h). This inequality means that 
the magnetization of each group of k spins is higher than 
that of k independent spins. Since the magnetization of 
different groups is additive, we conclude that the con- 
straint of tying together a group of spins increases the 
magnetization; therefore, the minimal M n (T,h) occurs 
at tt being the identity. 

Finally, the fact that the h = susceptibility exceeds 
the paramagnetic value follows from the two facts that 
M(T, 0) = and M(T, h) is greater than the paramag- 
netic value for all h. Q.E.D. 

Some Generalizations. — 

(A) The previous results hold true even if there are 
several species of bosons in the system, with possibly 
different masses, as long as there are no fermions present. 
The only difference is that the degeneracy of the ground 
state is the product of the degeneracies found in (i) , (ii) , 
(iii) above for each species. 

(B) The results hold for a 'relativistic' boson sys- 
tem obtained by replacing p 2 /2m by \/p 2 + m 2 = 



V^V 2 
of exp 



The underlying reason is that the kernel 



-/3V* 5 



Since K(X,X';T) > and > on a set of positive X 
measure, C w (T,h) > 0, and we have that M(T,h) > 



is also positive [[12], Eq. (7.11.11)]. 
Lest the reader think that positivity is a triviality, we 
remark that it does not hold for exp [— /3p 4 ] . 

(C) Our results hold in any dimension, of course. They 
also hold in lattice models (e.g., sometimes called 'Bose- 
Hubbard' models) in which the A in the kinetic energy, 
—A = —V 2 , is replaced by the discrete lattice Laplacian 
(= second difference operator) or Y^u j\ a \ a j m t ne usual 
second quantized notation. No regularity of the lattice is 
needed; indeed, the 'lattice' can be any general connected 
graph. 

Several recent works have discussed boson models with 
an antiferromagnetic coupling in the effective low-energy 
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Hamiltonian, resulting in a polar ground state. Our the- 
orem eliminates the possibility of the Hamiltonian (||) 
producing an antiferromagnetic state for bosons. Fur- 
thermore, it follows that the low energy parameteriza- 
tion (Q) of the Hamiltonian (|2|) must have C2 = 0. To 
be specific, we discuss the spin-1 case, but the situation 
is similar for bosons with higher spin. The potential (jl]) 
with C2 = satisfies the conditions of our theorem, and 
therefore its ground state for spin-1 bosons is degenerate, 
having states with J — N, N—2, N— 4, .... Switching on 
the spin-spin interaction, with non-zero C2, would choose 
between these degenerate states, break the degeneracy, 
and contradict the theorem. Therefore, a non-zero C2 
is inconsistent with a correct low energy theory of the 
Hamiltonian (0). Possible sources for such spin-spin in- 
teraction are then the (rather weak) direct dipole-dipole 
magnetic interaction, or exchange terms coming from the 
electron transfer between the atoms, thus going beyond 
the bosonic description of the system. 

In a recent experiment j?J the magnetic structure of a 
Bose-Einstein condensate of spin-1 87 Rb atoms was ex- 
plored. The authors studied the distribution of the dif- 
ferent magnetic number m values, but did not look at 
the global polarization properties. Based on the above, 
we suggest looking at the magnetization at finite field 
strength. According to the paramagnetic lower bound 
, the magnetization at the experimental temperature 
should exceed half of the maximal magnetic moment at 
fields as low as 1G. One can also rotate the macroscopic 
magnetic moment of the condensate to a specific direc- 
tion (by an appropriate magnetic field), and then (after 
switching off the field and allowing for relaxation) mea- 
sure the distribution of the m values, quantized along 
the same axis. A macroscopic moment should show up 
as a large fraction of the bosons in the m = 1 state. As 
discussed above, violation of the above predictions would 
indicate the importance of additional interaction terms, 
going beyond the Hamiltonian (||) . 

In conclusion, we have shown that fully polarized states 
are among the ground states of interacting bosons. More- 
over, the finite temperature magnetization and zero- 
field susceptibility were shown to be bigger than that of 
Maxwell-Boltzmann (distinguishable) free spins. Our re- 
sults hold in the presence of random on-site disorder and 
density-density interactions. They hold even if there is no 
interaction potential at all, because our theorems apply 
to the Hamiltonian H = ^p 2 /2m as well. They might 
be relevant for recent experiments suggesting the forma- 
tion of an internal magnetic moment at the onset of spin- 
triplet superconductivity, and to recent magnetic mea- 
surements of Bose-Einstein condensates. In particular, 
our ground state result constrains effective low-energy 



theories for interacting bosons, namely, it excludes the 
possibility of a spin-spin two-body term in the absence 
of explicit spin-dependent interactions in the underlying 
Hamiltonian. 

We thank Karsten Held, Robert Seiringer, and Jakob 
Yngvason for many useful comments. EHL thanks the 
National Science Foundation, grant PHY 0139984, for 
partial support of this work. 

Note added in proof.- After this paper was submitted 
for publication we learned that similar mathematical re- 
sults were obtained by Andras Siito Jl5| in a study of 
"cycle percolation" in Bose gases. Our theorem that the 
pure paramagnetic value is a lower bound to the magne- 
tization for nonzero field was not given there, however. 
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